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Abstract 

The scheme of Clauser and Dowling (Phys. Rev. A 53, 4587 (1996)) for factoring 
N by means of an A-slit interference experiment is translated into an experiment 
I , with a single Mach-Zehnder interferometer. With dispersive phase shifters the ra- 

(3J[)! tio of coherence length to wavelength limits the numbers that can be factored. A 

conservative estimate permits N ~ 10 7 . It is furthermore shown, that sine and 
cosine Fourier coefficients of a real periodic function can be obtained with such an 
Q\ ' interferometer. 
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Recently Clauser and Dowling (CD) have shown that factors of an integer N can be 
determined by simply measuring the peaks of the intensity distribution on the screen 
behind a Young's A-slit arrangement [|IJ . This device is distinct from the currently much 
debated proposals for quantum computation, because it does not need the entanglement 
of several quantal degrees of freedom. Therefore it is very immune against decoherence 
and could readily be implemented with present technology. The drawback is that it will 
not exhibit the potential increase in computational power characteristic of entanglement. 
Nevertheless, the work of CD indicates that single particle interference arrangements have 
useful applications beyond physical measurements. 

The purpose of this note is to point out that the CD-proposal can be translated into 
an experiment with only a single Mach-Zehnder interferometer. This will enhance the 
flexibility of this calculational device. A further point is that a Mach-Zehnder interferomter 
can also perform other computations, in particular cosine and sine Fourier transformations. 

Let us first focus on the work of CD. This proposal shows, that in a suitably chosen 
central region on the detection screen behind an A-slit arrangement, all intensity peaks are 
equal if, and only if, the quantity n = XR/a 2 is a factor of N, where A is the wavelength 
of the incident radiation, a is the center-to-center distance of the slits, and R is the 
distance between the slits and the screen. Furthermore, the Fraunhofer limit is assumed 
(R 3> a 3> A). Different values of n can be tested by adjusting any of the parameters, for 
instance A. The probability amplitude caused by slit i at a point r on the screen is given 
by ipi(r). The probability that the particle will hit this point is therefore 
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For the idea to be presented here, it is important to notice that the essential properties of 
I(r) can also be obtained in measurements with a series of different two-slit experiments. 
This is so, because quantum mechanical probabilities are the second power of probability 
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amplitudes. Therefore, the interference term of I(r), which is the double summation on 
the right hand side, only contains products made up of just two factors. And this feature is 
already present in the two-slit experiment. (This would not be the case, if, say, probabilities 
were the third power of probability amplitudes. Then the paradigm of quantum physics 
would only be fully contained in a three slit experiment.) The analytic behavior of the 
interference term of the iV-slit experiment can therefore be recaptured, when adding the 
intensities of all permutations of two-slit experiments obtainable by covering iV — 2 of the 
TV slits. This procedure will not reproduce the behavior of the mean intensity. But we 
can neglect the mean intensity here, because it is only of statistical relevance. 

Instead of doing several two-slit experiments with different distances of the slits, it is 
of course more convenient to use a Mach-Zehnder interferometer and vary the phase shift. 
A further advantage of a Mach-Zehnder interferometer is that it has only two outputs, 
instead of the continuum of a screen. Then, with an arrangement as shown in Fig.l, 
factoring can be achieved in the following way: The phase shift x is increased in discrete 
steps 2n/n, such that Xj — 2irj /n. For simplicity we assume that at each step one particle 
is sent into the interferometer. We look at the outputs only when j = kN, k = 1, 2, 3... . 
So at the k th observation the probability to register the particle at detector A is 



We have Pa{^) = 1 only if n is a factor of N. Suppose we perform n observations, starting 
with k — 1 until k = n. Then the sum of the particles registered in detector A will be 



For n a factor of N we have I n = n, because then the phases in the cosine terms will all 
be multiples of 2n. However, if n is not a factor, the cosine-terms of the summation will 
roughly cancel each other, and we will only have I n ~ n/2. To see this more clearly, we 
write N/n = M + r, where M is the largest integer for which N/n > M, and r is the 
rational remainder, r = L/n, with L e [l,2,....,n — 1]. Then the summation over the 
cosine terms can be written as 



The phases of the cosines are now equidistantly spread in n steps from the value 2irL/n, 
which is less than 2tt, to 2nL, which is at most 27r(n — 1), so that when reducing the 
phases to the interval [0, 2n], this interval will be used quite evenly. 

For practical applications it is important to know what calculations are possible. Since 
the best interferometers today are operated with light, it is sufficient to look at the char- 
acteristics of such interferometers. If we implement the phase shift in the usual manner 
by increasing the path length difference of the two arms of the interferometer, the coher- 
ence length of the incident radiation sets the limit to the largest number N that can be 
factored. According to eq.(3) the maximum phase shift to be set in this scheme is 2irN. 
Assuming a wavelength of A = 500nm and conservatively limiting the coherence length, 
given by C = A 2 /AA, where AA is the standard deviation of the wavelength distribution, 
to C — 5m, numbers up to iV 10 7 can be factored. 
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One also wants to know how long a calculation will take. When testing whether n is a 
factor of N, one must set nN different phase shifts. Since the largest factor to be checked 
is of the order of y/~N, the longest check will take a time proportional to N 3 ^ 2 . 

A further consideration is the maximum deviation permissible in the phase increments 
2n/n, in order to ensure correct identification of the factors of N. Suppose the actual 
increment is 2ir/(n + d), where d is the deviation. In the ideal case we have d — 0, and 
when n is a factor of N the sum (3) yields an intensity of n. Now this sum will still be 
substantially larger than n/2 (which would indicate a non-factor), if the deviation d is 
such, that for the largest phase shift to be set (at k = n) we have 2nN — | rather than the 
ideal value 2irN. Then the contribution of the last term of the sum is zero, whereas all the 
others are still positive, so that the total intensity will be roughly n(| + -). This limits 
;he permissible relative deviation of the phase shift increments from their ideal value to 
- < (AN) -1 . For our example of above this would mean - < 2.5 x 10 -8 . Such accuracy 
can be achieved in optical interferometers, when the lengthening of one arm relative to 
the other is itself controlled interferometrically. 

It is also possible to implement a kind of parallel computation with a setup as illustrated 
in Fig. 2. Here, the numbers rii, n 2 , n-j can almost simultaneously be checked for being 
factors of N, thereby utilizing an incident particle for more than just one computation. 
The phase shifters in the various loops are simultaneously stepped up, but at the different 
increments 2n/rii, 2n/n 2 , 2n/n 7 , respectively. Let us assume < rij for i < j. And 
for the sake of convenience we sum the intensity in a given detector as needed for testing 
for the largest of the possible factors of the loops involved. For instance, at detector A 
(and also at B) we would sum until k = n 4 , thus 
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The expectation values for I a, Ib and Ic+Id (the latter also just summed up to 114) for the 
various possibilities are shown in Table I. Clearly, the results permit unique identification 
of the eight possible answers of interest. 



Table I. Dependance of intensities on property of "factor" or "non-factor". 
Factors are indicated by "F", non-factors by "— ". Intensities in units of 714/8. 
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The situation is a little bit more complicated for the interferometer loops following the 
lower arm of the ni-interferometer. We must account for the possibility that n\ is a factor 
of N. If we then measured the intensity at any of the detectors E to H after every iV 
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phase increments, we would never detect a particle. Therefore it is necessary to start the 
loops on the lower arm of n\ with a delay of around ni/2 increments. Similar arguments 
apply to any of the interferometers connected to an output which is dark when the phase 
shift is a multiple of 2tt. For instance, testing whether n.3 and are factors of N requires 
that one waits for about (rii + n$)/2 increments, and from then on sums the intensities at 
E (or at F) after every N further increments. Then the expectation values for Ie for the 
eight possibilities (i.e., which of n 1 , n 3 , n Q is a factor of N) are analogous to those shown 
for I a in Table I. However, the offset relative to the total number of phase increments to 
be gone through is small, such that one can still speak of essentially parallel computation. 
But, of course, as this is a device not exploiting the entanglement of quantum systems, the 
number of loops, or the total amount of time, needed for factoring N rises polynomially 
with N, rather than just logarithmically. 

Let us now turn to how Fourier transformations can be performed by means of an 
interferometer. Specifically, it is possible to obtain the cosine and sine Fourier coefficients 
of a real valued positive function /(£), which is periodic with the period r, and where 
t denotes time. For this it is necessary that the intensity incident on the interferometer 
varies in time proportional to f(t), and that the phase shift \ applied between the two 
paths of the interferometer increases linearly with time, e.g. 

„ , 2irmt , , 

Xif) = — , (6) 

where m is an integer. In order to avoid problems of the change of energy of the particles 
when experiencing a time-dependent phase shift, we limit ourselves to slow rates, m/r <C 
E/h, where E is the energy of the particles and h is Planck's constant. Then the difference 
of the intensities at detectors A and B integrated over one period r is 
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cf f{t) cos (^)dt, (7) 



where c is the proportionality constant between f(t) and the incident intensity. As can 
be seen, I a — Ib is proportional to the m th cosine Fourier coefficient of /(£). In order to 
obtain the sine Fourier coefficients, one only needs to add a constant phase shift of —n/2 
to xif). 

Note that this calculational procedure is different from the experimental method em- 
ployed in Fourier spectroscopy, which has been used in astronomy since the days of Michel- 
son and which has recently also been used in electron B and neutron || interferometry. In 
Fourier spectroscopy one can measure the spectral distribution of the incident radiation, 
if this distribution is symmetric about the mean, but one gets no phase information. 
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Figure Captions 

Fig.l: Ideal Mach-Zehnder interferometer with semitransparent mirrors for beam 
splitting at the entrance and beam superposition at the exit. A phase shifter introduces 
the additional phase x between the two paths. Particles are registered in detectors A 
and B. 

Fig. 2: Series connection of several Mach-Zehnder interferometers for almost simulta- 
neous determination whether the integers rii, ... , are factors of N. The particles 
are registered at detectors A to H . 
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